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INTRODUCTION 


This report, prepared under Contract NAS2-10331 (MOD 10), outlines the devel- 
opment and validation of an interface between the three-dimensional transonic 
analysis program R0T22 (Reference 1) and the DATAMAP (Reference 2) progreun. 

Following the section describing the development of the interface, the valida- 
tion is carried out as follows. First, the DATAMAP program is used to analyze 
the instrumented AH-lG flight test data (Reference 3). Specifically, records 
2872 and 2873 are analyzed at an azimuth of 90**, and Record 2806 is analyzed 
at 60^. Trim conditions for these flight conditions are then calculated using 
the Bell performance prediction program ARAM45 (Reference 4). Equivalent 
shaft, pitch, and twist angles are calculated from ARAH45 results emd used as 
input to the R0T22 program. The interface uses the R0T22 results and creates 
DATAMAP information files from which the surface pressure contours and sec- 
tional pressure coefficients are plotted. Twist angles input to R0T22 program 
are then iteratively modified in the tip region until the computed pressure 
coefficients closely match the measurements. 

The results of this match are discussed. Conclusions derived from this work 
are then summarized. 
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DEVELOPMENT OF THE INTERFACE 


ROT22 was modified to pass its results to DATAMAP. The appropriate R0T22 re- 
sults at different station and chord values are interpolated to a rectangular 
grid usable to DATAMAP. The interpolated values are then written to a file in 
a format appropriate for DATAMAP to process. Variable identification, plot 
identification, plot titles, and the grid definition are also written to the 
file. DATAMAP then processes the file and stores the information for further 
analysis and plotting. 

Program ARSA05 is the result of the above development. A user's guide to this 
program is contained in the appendix which details the program's function, in- 
put information, and output description. 




3 



OUTLINE OF THE INTERFACE 


ROT22 has been modified to write Data Transfer Files (DTFs) in the format 
specified in USAAVRADCOM-TR-80-D-30B. Now that a connection to DATAMAP 
through a DTF is established, DATAMAP is able to produce the required plots. 
However, several analytic and procedural abilities (in addition to writing 
DTFs) were added to ROT22 to produce the data required by DATAMAP. These in- 
clude (see Figures lA through 1C): 

a. Conversion of static pressure coefficient to local absolute pres- 
sure. 

b. Chordwise and spanwise interpolation of blade thiclcness, local Mach 
numbers and local absolute pressure. 

c. Conversion of local absolute pressure to static pressure coeffi- 
cient. 

d. Data extension to obtain pseudo rotor cycles. 

First, DATAMAP cannot compute backwards from the static pressure coefficient 
(Cp) to cdssolute pressure so this analysis was added to R0T22. Thus, R0T22 
writes negative Cp, local absolute pressure and incidentally local Mach number 
on the DTF. Second, R0T22 computes values at station (z direction) and chord 
(y direction) values that are, in general, different from the values desired 
in DATAMAP. R0T22 will not provide values at a user specified grid via two- 
way interpolation. Third, the local Cp is computed from the interpolated 
local absolute pressure since Cp cannot be meaningfully interpolated. Fourth, 
R0T22 computes values for one azimuth during a run. To enable DATAMAP to use 
azimuth for computations and displays, R0T22 provides constant-valued time 
histories for each of the three quantities at each chord (y) and span (z) 
position that includes one complete rotor cycle. The time histories will be 
correct only at the specified azimuth for which R0T22. computations were ac- 
tually performed. However, the simulation of a complete rotor cycle allows 
DATAMAP to use azimuth in computations and displays of the data. 


5 


pqccCDMQ PAGE 91ANK WT FBJSP 




6 







7 











ROT22 

RESULTS 


Find cross section 
for intermediate 
interpolation 


Compute local 
absolute pressures 


b. 


Interpolate blade 
thickness at the 

same 

points for 

each 

cross sec- 

tion 

( rectang- 

ular 

grid) 


Interpolate local 
Mach numbers at 
the same points 
for each cross 
section (rectang- 
ular grid) 


Interpolate local 
absolute pressure 
at the same points 
for each cross 
section (rectang- 
ul 


Interpolate blade 
thickness from 
rectangular grid 
to user specified 
grid (2-way inter- 
polation) 


Interpolate local 
Mach numbers from 
rectangular grid 
to user specified 
grid (2-way inter- 
olatlon) 


Interpolate local 
absolute pressure 
from rectangular 
grid to user spe- 
cified grid (2-way 
interpolation) 


Write azimuth and 
Mach number item 
code records and 
data records on 
DTF 


Write local absolute 
pressure Item code 
records and data 
records on DTF 


Compute negative CPs 
at user specified 
rid 


1 d. 



Write local 

nega- 

tive CP 

item code 

records 

and 

data 

records 

on 

DTF 


Figure 1C. Information flow chart - R0T22 results 
(see text for explanation of a-d) 
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R0T22 writes five record types to the DTP. Three of the five record types are 
written in the input routines. These are Instruction Records, a Counter Re- 
cord, and Info File Records. After R0T22 is otherwise finished, the other two 
record types are written. These are Item Code Records and Data Records. 

The R0T22 results to be transferred to DATAMAP are the values on the blade 
surface. The quantities of interest are the station (z) and chord (x) values 
used in the R0T22 analysis, the station Mach numbers, local Mach numbers, 
local CpS and blade surface (y) at each point (measured from reference line). 
The CpS are converted to local absolute pressures by the relation 

local absolute pressure = [1+0. 7*Cp*( station Mach number)^]* 

ambient air pressure 

A rectangular grid is generated by chordwise interpolation of the surface (y) , 
local Mach number and local cdjsolute pressures. These values are then con- 
verted to the user specified grid via a two-way interpolation on the station 
and chord values. If, during the above interpolations, a point is extrapo- 
lated, a warning message is written. The local CpS are computed from the 
local absolute pressures by the relation 

Cp = [(local absolute pressure/ambient pressure) -l]/[.7* 

(station Mach number)^] 

The negative of Cp is passed to DATAMAP. 

After all the interpolations are finished, the Item Code Records and Data 
Records are written to the DTP. The azimuth and Mach numbers are written to 
the DTP. Five of these Data Records are written with the only difference 
being the azimuth position (0®, 90®, 180°, 270®, 360°). This generates the 
pseudo time history. Next the local absolute pressures are written to the DTP 
five times. Finally, the local negative CpS are written to the DTP five 
times. 
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THEORETICAL CONSIDERATIONS 


Various equations used to calculate pressure coefficients are provided in this 
section. In addition, equations used for calculating the equivalent shaft 
angle, pitch angle, and twist angles are also developed here. 


Pressure Relationships 


C 


P 



( 1 ) 


where is coefficient of pressure, p^ is absolute pressure, p^^^ is ambient 
pressure, p is air density, is sectional free stream velocity defined as 

= Qr + Vp sin «1> (2) 

where Q is angular velocity of the blade, r is the radial location, is 
flight velocity, and <ii is blade azimuth. 


From gasdynamic relations equation (1) can also be written as 


C 


P 


ym: 


Pop 


1 


(3) 


or 


C 


P 






1 


(4) 


where is the sectional free-stream Mach number, y is the ratio of specific 
heats, a^ is speed of sound in the ambient, and a is the local speed of sound 
in the flow. 
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In DATAMAP program, Eq. 1 is used where the ambient pressure and density are 
input. On the other hand, Eq. 4 is used in ROT22 program where the ambient 
speed of sound is input. In any event. Equations 1, 3, and 4 are identical 
when perfect gas assumptions are made. 

Flight Condition Simulations 


Changes in the blade sectional angles of attack due to flapping, pitch- flap 
coupling, and blade cyclic are not represented in R0T22 analysis. However, by 
equating the relationship for the effective angle of attack in R0T22 with that 
in a performance analysis such as ARAM45 program, one can readily obtain rela- 
tionships for equivalent twist and pitch £mgles. 

Similarly, one can obtain an equivalent shaft tilt angle. In R0T22 program. 


In ARAM45 progreun. 


I 



c 


r 


a 


TPP 


a - a 


1 


( 6 ) 


where ot^pp is tip path plane tilt with respect to flight velocity vector, 
is the shaft tilt angle, and a^^ is Fourier coefficient from blade flapping 
angle written as 


P = a^ + aj^ cos i)> + bj^ sin ^2 


(7) 


In the present analysis, shaft tilt aft is defined to be positive. 


Upon equating equations 5 and 6, 


^“s^Eq. ■ ^“s ■ ®1^ARAM45 


( 8 ) 
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Both and on the right-hand side of Eq. 8 can be obtained from ARAM45 
results . 


Equivalent pitch and twist: 

Blade element angle of attack relationship in ROT22 can be written as. 


a = 0 + 

r o 


e„ + tan 


-1 


AflR 


Qr + 


Vp sin »ji 


(9) 


The same relationship in ARAM45 analysis is, 

o(^ = 6^ + 6^ - Bj^ sin i(i - cos »|» - P tan 6^ + <|> (10) 

where 6^ is collective pitch, 6^ ^ is geometric twist at radius r, Bj^, and A^^ 
are cyclic pitch coefficients, 6^ is pitch flap coupling, A is inflow ratio, 
and ^ is inflow angle defined as 


tan <)i = 


•r p sin - Vp cos cos 4> sin p + (Vp sin - v^) cos p 

Qr + Vp sin «1» 


( 11 ) 


where p is first azimuthal derivative of blade flapping motion and is in- 
duced velocity. Since ARAH45 program calculates <|>, one need not use Eq. 11. 
On the other hand, R0T22 program does not e:q>licitly solve for the inflow 
ratio. Instead, it solves finite difference flow field equations to calculate 
local velocity components. Though one can use these to calculate an effective 
inflow ratio, it is a very time consuming and computationally demanding task. 
Here it is assumed that 


^ARAM45 ” \oT22 

Now Eq. 9 can be equated with Eq. 10 to get 

6 = 6 - B, sin >if - A, cos »1> - P tan 6, 

o_ o 1 ^ 1 -r f 

Eq 


( 12 ) 


(13) 
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( 14 ) 


6m = 6m + 4> ■ tan ^ — — 

T T ^ Qr + 

Eq 


where 

t|) = azimuth angle 

0 = the equivalent pitch angle 

°Eq 

0^ = the equivalent twist at station 


V_ sin «1» 
r 



ANALYSIS 


In this report, the R0T22 - DATAMAP interface is validated by using three 
level flight records acquired by NASA-ARC of the instrumented AH-lG helicop- 
ter. This instrumented helicopter has a two-bladed main rotor with rectangu- 
lar planform blades of 28.63 inch chord. The airfoil section used on this 
blade is a symmetrical section of 9.3% thickness. Airfoil coordinates for 
this blade are provided in Table 1. 

Flight conditions representing records 2806, 2873, and 2872 are provided in 
Table 2. The ARAM45 results, equivalent shaft angle, equivalent pitch emgle, 
and equivalent twist angles, are provided in Table 3. 

These equivalent angles are input to R0T22 program to calculate blade surface 
pressure distribution for each of the records selected. 

The normal force coefficients computed for these records are compared in Fig- 
ures 2 to 4 with the normal force coefficients derived from flight tests. In 
these figures, the flight test data is shown by a solid line and the R0T22 
data by a broken line. Initial plans called for approximately matching these 
normal force coefficients by modifying the tip region twist. However, a 
careful examination has revealed that the combination of test data scatter and 
the linear interpolation/extrapolation scheme used in DATAMAP program makes 
the normal force calculations unreliable. Consequently, tip region twist is 
modified to approximately match the pressure coefficients rather than the 
normal force coefficients. It will be shown later that the agreement between 
the normal coefficients worsens when the agreement between the pressure coef- 
ficients is improved. Since we already discussed that the normal coefficients 
calculated from test data are unreliable, we may not have to be concerned with 
this discrepancy. 
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TABLE 1. AIRFOIL COORDINATES 


0.0 

0.0 

0.00453 

0.01163 

0.00862 

0.01569 

0.01250 

0.01857 

0.01628 

0.02087 

0.02366 

0.02453 

0.03092 

0.02743 

0.03810 

0.02986 

0.04524 

0.03195 

0.05590 

0.03463 

0.06652 

0.03689 

0.08064 

0.03939 

0.09121 

0.04098 

0.10177 

0.04235 

0.11231 

0.04353 

0.12285 

0.04455 

0.13339 

0.04543 

0.15093 

0.04660 

0.17548 

0.04773 

0.20001 

0.04836 

0.21402 

0.04851 

0.22102 

0.04855 

0.22452 

0.04855 

0.23503 

0.04851 

0.25253 

0.04832 

0.28053 

0.04767 

0.30152 

0.04696 

0.32950 

0.04572 

0.35048 

0.04461 

0.38544 

0.04244 

0.42036 

0.04015 

0.45528 

0.03785 

0.49020 

0.03556 

0.52512 

0.03326 

0.56003 

0.03097 

0.59495 

0.02867 

0.62987 

0.02638 

0.66479 

0.02408 

0.69970 

0.02179 

0.73462 

0.01949 

0.76954 

0.01720 

0.80446 

0.01490 

0.83938 

0.01261 

0.87430 

0.01031 

0.90921 

0.00801 

0.94413 

0.00572 

0.97905 

0.00342 

1.00000 

0.00205 




TABLE 

2. FLIGHT CONDITIONS 


1 — 1 

Rec # 

2872 

2873 

2806 

V_, KTAS 
r 

150.2 

143.7 

123.7 

N, rpm 

301 

315.5 

316 


90® 

90® 

o^ 

o 

o 

SW 

8150 

8050 

8100 

OAT °F 

68.9 

63.5 

69.1 

Pq* Psi 

11.98 

11.56 

13.17 

0 

0 

18.594° 

16.796® 

14.92 

h' 

0°, 6.405® 

0®, 4.8839® 

0, 2.808 

a 

o 

2.75® 

2.75® 

2.75 

“r “i 

-2.385®, -1.61® 

-1.6108®, -1.1737® 

-1.743, -0.975 


o 

o 

o 

o 

O 

o 

o 

o 

o 

o 


0® 

0® 

0® 

Of 

s ' 

-6.7144® 

-6.365 

-5.87 

1 

-0.052378® 

-0.0467 

-0.03902 

1 ^ — - - 1 
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TABLE 3. CALCULATED EQUIVALENT CONDITIONS 


Rec # 

2872 

2873 

2806 

a 

hq 


-4. 

3294 


-4.755 

- 

■4.13 

0 

°Eq 


12 

.185° 


11 

.91° 

12.488 

r/R = 

T 

r 

<t>“ 

(e„ )« 

r Eq 

u 

H 

CD 

(|,0 

(e„ ) 

r Eq 

r 

(0T ) 

r Eq 

0.15 

-1.5 

-6.6 

-2.25 

-1.5 

-6.1 

-2.05 

-6.5 

-2.36 

0.3 

-3.0 

-5.65 

-4.11 

-3 

-5 

-3.75 

-4.95 

-3.85 


-5.0 

-4.9 

-6.4 

. 

-t) 

-4.2 

-5.97 

-3.9 

-5.89 

0.75 

-7.5 

-4.35 

-9.14 

-7.5 

-3.62 

-8.63 

-3.2 

-8.45 

0.85 

-8.5 

-4.2 

-10.2 

-8.5 

-3.5 

-9.7 

-3.03 

-9.48 

0.9 

-9.0 

-4.13 

-10.74 

-9.0 

-3.4 

-10.2 

-2.95 

-9.99 

0.95 

-9.5 

-4.06 

-11.26 

-9.5 

-3.325 

-10.73 

-2.86 

-10.48 

1.0 

-10.0 

-4.0 

-11.78 

-10.0 

-3.2674 

-11.25 

-2.79 

-10.99 
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Figure 2. Normal force coefficient for Record 2872. 
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Figure 3. Normal force coefficient for Record 2873. 
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Figure 4. Normal force coefficient for Record 2806. 
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RESULTS 


In this section, the sectional pressure coefficients calculated by R0T22 are 
compared with those from flight tests. Then the sectional pressure coeffi- 
cients and surface pressure contours obtained after modifying the tip region 
twist to approximately match the pressure coefficients are provided for each 
of the records. In addition, measured surface pressure contours are also pro- 
vided. Finally some discussion of these results is provided. 

Record 2872 


Sectional pressure coefficients calculated by ROT22 program using the equiva- 
lent angles provided in Tcd>le 3 are compared with the flight test results in 
Figures 5 through 9. From these it was concluded that the upper surface suc- 
tion was underestimated considerably. This mandated an increase in the twist 
in the tip region. After a few iterations, the new twist rate provided in 
Table 4 was determined. The results from R0T22 after modifying the tip region 
twist angles are shown in Figures 10 through 16. Here, the pressure contours 
from ROT22 results shown in Figures 15 and 16 may be compared with those from 
flight tests (Figures 17 and 18). The normal force coefficients calculated 
after modifying the twist are provided in Figure 19. As mentioned in the 
earlier section, the agreement between the theory and flight tests has wors- 
ened when compared with Figure 2. 


Record 2873 


Results for Record 2873 are provided in a manner identical to Record 2872. 
Here the sectional coefficients of pressure before adjusting the twist are 
provided in Figures 20 through 24. The same after adjusting the twist are 
shown in Figures 25 through 29. 

The surface pressure contours from R0T22 results are provided in Figures 30 
and 31. These may be compared with those in Figures 32 and 33 where the 
flight test results are used. Finally, the normal force coefficients obtained 
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Figure 5. Sectional pressure coefficient at 86% radius for Record 2872. 
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Figure 6. Sectional pressure coefficient at 91% radius for Record 2872. 
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Figure 7. Sectional pressure coefficient at 95% radius for Record 2872. 
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Figure 8. Sectional pressure coefficient at 97% radius for Record 2872. 
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Figure 9. Sectional pressure coefficient at 99% radius for Record 2872. 
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TABLE 4. TWIST RATES, 

(0.T. )r.rt/ TO MATCH SURFACE 

' Tj^'EQ' 

PRESSURES 


r/R Record 2872 

Record 2873 

Record 2806 


0.75 

-9.37 

-8.63 

-8.45 

0.85 

-10.20 

-10.0 

-9.50 

0.90 

-10.74 

-10.40 

-9.70 
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-10.50 

-10.10 

• -9.10 
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-11.00 

-10.0 
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Figure 10. Sectional pressure coefficient at 86% radius 
for Record 2872 (modified twist) . 
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Figure 11. Sectional pressure coefficient at 91% radius 
for Record 2872 (modified twist). 
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Figure 12. Sectional pressure coefficient at 95% radius 
for Record 2872 (modified twist). 
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Figure 13. Sectional pressure coefficient at 97% radius 
for Record 2872 (modified twist) . 
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Figure 14. Sectional pressure coefficient at 99% radius 
for Record 2872 (modified twist) . 
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Figure 15. R0T22 blade tip pressure contours for 

Record 2872 (upper surface). 
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Figure 16. R0T22 blade tip pressure contours for 

Record 2872 (lower surface). 
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Figure 17. Flight test blade tip pressure contours 
for Record 2872 (upper surface). 
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Figure 18. Flight test blade tip pressure contours 
for Record 2872 (lower surface). 
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Figure 19. Normal force coefficient for Record 2872 
(modified twist) . 
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Figure 20. Sectional pressure coefficient at 86% radius for Record 2873. 
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Figure 21. Sectional pressure coefficient at 91% radius for Record 2873. 
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Figure 22. Sectional pressure coefficient at 95% radius for Record 2873. 
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Figure 23. Sectional pressure coefficient at 97% radius for Record 2873. 


43 









NEG CP 



FRPICTN OF CHORD 


o o o 

DERIVED PARAMETER: 

2873 GROSS WT 

R/RADIU5 LONG CG 

BLADE STATIC PRESSURE 

COEFF 

SHIP MODEL AH-IG 
UPPER SURFACE 

A 

A ^ COUNTER 

2873 

GROSS WT 


SHIP MODEL 

AH-IG 


.99 

R/RADIUS 

LONG CG 


LOWER SURFACE 

1 DERIVED PARAMETER: 

BLADE STATIC PRESSURE 

COEFF 




COUNTER 

441053 

GROSS WT 

8050 

SHIP MODEL 

OLSAHl 


.99 

R/RADIUS 

LONG CG 


UPPER SURFACE i 

CYCLE 

AVERAGE; 

LOCAL NEG CP: RQT22 





CnUNTER 

441053 

GROSS WT 

8050 

SHIP MODEL 

OLSAHl 


.99 

R/RADIUS 

LONG CG 


LOWER SURFACE 

CYCLE 

AVERAGE: 

LOCAL NEG CP: RDT22 



i 


BHT,U5ftRTL DATAMAP (VERS 3.10 - 24/02/82) 04/16/83 BELL HELICOPTER 


Figure 24. Sectional pressure coefficient at 99% radius for Record 2873. 
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Figure 25. Sectional pressure coefficient at 86% radius 
for Record 2873 (modified twist). 
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Figure 27. Sectional pressure coefficient at 95% radius 
for Record 2873 (modified twist). 
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Figure 28. Sectional pressure coefficient at 97% radius 
for Record 2873 (modified twist). 
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Figure 29. Sectional pressure coefficient at 99% radius 
for Record 2873 (modified twist). 
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Figure 30. 


R0T22 blade tip pressure contours for 
Record 2873 (upper surface). 
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Figure 31. R0T22 blade tip pressure contours for 

Record 2873 (lower surface). 
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Figure 32. Flight test blade tip pressure contours 
for Record 2873 (upper surface). 
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Figure 33. Flight test blade tip pressure contours 
for Record 2873 (lower surface). 
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after modifying the tip region twist input to ROT22 are compared with those 
from test data in Figure 34. As in the case of Record 2872, the agreement 
becomes worse when compared to the results shown in Figure 3. 

Record 2806 


The results here are presented for an azimuth angle of 60® in a manner iden- 
tical to the preceding records. In other words, sectional pressure coeffi- 
cients from ROT22 program are first provided in Figures 35 through 39. The 
twist input to R0T22 is then modified to improve pressure coefficient correla- 
tion in the tip region. The ensuing results are then provic.ed in Figures 40 
through 46. While Figures 40 through 44 provide the sectional pressure coef- 
ficients, Figures 45 and 46 provide computed surface pressure contours. These 
may be compared with the measured surface pressure contours provided in Fig- 
ures 47 and 48. Finally, the normal force coefficients calculated from ROT22 
results are compared with measurements in Figure 49. 
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Figure 34. Normal force coefficient for Record 2873 
(modified twist). 
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Figure 35. Sectional pressure coefficient at 86% radius for Record 2806. 
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Figure 36. Sectional pressure coefficient at 91% radius for Record 2806. 
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Figure 37. Sectional pressure coefficient at 95% radius for Record 2806. 
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Figure 38. Sectional pressure coefficient at 97% radius for Record 2808. 


59 




NEG CP 


O 



FRACTN OF CHORD 


O O O 

DERIVED PARAMETER; 

2806 GROSS WT 

R/RADIU5 LONG CG 

BLADE STATIC PRESSURE 

COEFF 

SHIP MODEL AH-IG 
UPPER SURFACE 


A A COUNTER 

2806 

GROSS WT 


SHIP MODEL 

AH-IG 


^ .99 

R/RADIUS 

LONG CG 


LOWER SURFACE 

DERIVED PARAMETER: 

BLADE STATIC PRESSURE 

COEFF 




COUNTER 

523565 

GROSS WT 8100 

SHIP MODEL 

OLSAHl 


.99 

R/RADIUS 

LONG CG 


UPPER SURFACE 

CYCLE 

AVERAGE: 

LOCAL NEG CP; R0T22 





COUNTER 

523565 

GROSS WT 

8100 

SHIP MODEL 

OLSAHl 


.99 

R/RADIUS 

LONG CG 


LOWER SURFACE 

CYCLE 

AVERAGE: 

LOCAL NEG CP; R0T22 





BHT.USARTL DfiTAMflP (VERS 3. 10 - 24/02/82 ) 04/16/83 BELL HELICOPTER 


Figure 39. Sectional pressure coefficient at 99% radius for Record 2806. 
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Figure 40. Sectional pressure coefficient at 86% radius 
for Record 2806 (modified twist). 
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Figure 41. Sectional pressure coefficient at 91% radius 
for Record 2806 (modified twist). 
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Figure 42. Sectional pressure coefficient at 95% radius 
for Record 2806 (modified twist) . 
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Figure 43. Sectional pressure coefficient at 97% radius 
for Record 2806 (modified twist). 
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Figure 44. Sectional pressure coefficient at 99% radius 
for Record 2806 (modified twist) . 
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Figure 45. 


R0T22 blade tip pressure contours for 
Record 2806 (upper surface). 
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Figure 46. R0T22 blade tip pressure contours for 

Record 2806 (lower surface). 
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Figure 47. Flight test blacie tip pressure contours 
for Record 2806 (upper surface). 
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Figure 48. Flight test blade tip pressure contours 
for Record 2806 (lower surface). 
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Figure 49. Normal force coefficient for Record 2806 
(modified twist). 
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CONCLUSIONS 


In addition to demonstrating the capability of the interface, the following 
conclusions are drawn from this study. 


1 . 


In all the cases studied, the location of the shock is well predicted. 
However, the negative pressure coefficients are underpredicted. This 
indicates an apparent thickening of the blade section. This disagreement 
appears to be a meuiimum at 86% radius and diminishes as the tip is ap- 
proached. Since the instr\imented AH-16 blades were built to close toler- 
ance, this problem is not related to blade construction. It has been 
shown in References 5 and 6 that such increases in the negative pressure 
coefficients occur when the blade encounters a vortex. Also, the two- 
dimensional, unsteady transonic results of Reference 5 show that the in- 
creased suction on airfoils cem last well after the vortex has swept past 


Figure 50 shows the planfcrm view of a rigid wake trajectory 
for two flight conditions. The first one closely represents the flight 
conditions of Records 2872 and 2873. The second case represents the 
flight conditions of Record 2806. From both the figures, it may be noted 
that the blade region has just emerged from a vortex encounter within the 
preceding 15^* azimuth. Consequently, the blade-vortex interaction may be 
the source of increased negative pressure coefficients on the blade. 
Further investigation is recommended before emy final conclusion can be 
made. 


2. Use of R0T22 program in concert with performance prediction programs pro- 
vides a means of obtaining realistic pressure distribution and the loca- 
tion of shock in the tip regions of the blade. 


3. The scatter in the measured pressure distribution, together with the lin- 
ear interpolation/ extrapolation scheme used in the DATAMAP program, makes 
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Figure 50. Planform view of tip vortex trajectories for two 
flight conditions using rigid wake assumptions. 
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the normal force coefficient calculations highly unreliable. This prob- 
lem can be alleviated by smoothing the measured pressure data (while 
retaining discontinuity due to shock) and also by providing an option to 
insert data points near the trailing edge of the blade when measurements 
are faulty in that region. 
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APPENDIX 


User's Guide to Program ARSA05 

Program ARSA05 is Bell's version of the original NASA-Ames R0T22 program. It 
is a 3-D transonic, full potential, nonconservative, quasi-steady, inviscid 
analysis program developed to predict the pressure and flow velocity on heli- 
copter rotor blades. This version includes an interface to DATAHAP, a data 
management and analysis package, to provide extensive plotting capabilities. 

FUNCTION 

ARSA05 can be used to assess the effect of compressibility on various blade 
configurations and its implication on the aero-acoustic properties of the 
blade. The program can also be used to locate the high conqpressibility areas 
near the blade tip where transonic airfoils can be designed to minimize com- 
pressibility using two-dimensional design techniques (ADAH system) . Using the 
DATAHAP^ interface, the user cam obtain surface Mach number and pressure con- 
tour plots. In addition, sectional pressure coefficients, pressure and Mach 
number plots can be obtained. 

The user should take note of the following: 

1. The analysis is quasi-steady and, hence, it does not account for the time 
history of the blade motion. 

2. The analysis does not account for the blade flapping velocities or dy- 
namic twist. However, the blade azimuth location, blade deflection due 
to flapping, geometric twist, cyclic and collective pitch, blade contour 
and tip path plane inclination, with respect to the flight path, can be 
accounted for. 
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3. The wake is not represented accurately in the program because of the com- 
plex transformations involved. Consequently, the induced velocities are 
in error. The error would probably be larger xmder hovering conditions. 
Some empirical technique may have to be used to get accurate inflow 
through the rotor. Though this is a certain limitation for performance 
prediction, it is not anticipated to be a serious handicap in assessing 
the effect of compressibility on different blades. This is particularly 
true for helicopters in high-speed flight for which the blade pitch 
angles are small in the advancing blade tip region. 
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INPUT INFORMATION 


FORMAT: 

CARD 0 : 
CARD 1 : 
CARD 2: 


All the odd-numbered cards are skipped. The format for all the 
even-numbered cards is 8F10.0 unless otherwise noted. 

Title Card (format 20A4) 

Comments 

FNX, FNY, FNZ, PLTCP, PLTSEC, PLTSRF, FCONT, FPRINT 

FNX No. of chordwise grids in the entire flow field (use 30; Max = 
120 ) 

FNY No. of y grids (+ up) in the entire flow field (use 4; Max = 
16 ) 

FNZ No. of z grids (+ sp£uiwise out) in the flow field (use 8; Max - 
32) 

PLTCP controls CP plot routine, currently unavailable (Use 0.). 

PLTSEC controls sectional characteristics computation for sectional 
plots (Use 1.). 

PLTSRF controls surface characteristics computation for surface 
plots (Use 1.). 

FCONT: This sets the switch for restart options. 

0: Initial velocity potential, ^ values are generated and 

final 0 values are not saved (i.e.. Tape 4 and Tape 7 are 

not used) . 

1: Initial values are generated (Tape 4 en^ty) and final 0 

values are saved on Tape 7 for future use. 

2: Initial values are read in from Tape 4 from a previous 

run and the final 0 values are written out on Tape 7. 

3: Initial values are read in from Tape 4 but the final ^ 

values are not saved (Tape 7 enqpty). 

FPRINT controls printout of intermediate mesh size and final mesh 
size results. 

0: All sectional printouts suppressed for intermediate mesh 

sizes. 

?t0:. All sectional data printed. 
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CARD 3: 


Comments 


CARD 4: 


CARD 4a ; 

CARD 5 ; 
CARD 6: 


FITMAX, FITMIN, CONV, PlO, P20, P30, BETAO, FHALF 

FITMAX: Max no. of iterations in the first grid defined by CARD 2. 

(Use 10 for rough grid, 20 for intermediate, and 30 for 
fine mesh. ) 

FITMIN: Min no. of iterations in the first grid (Use 1.). 

CONV: Convergence criteria; (suggested 0.00003) 

PIO Subsonic relaxation parameter 1.0<P10<2.0 (use 1.75) 

P20 Subsonic relaxation parameter 0<P20<1 (Use 1.0) 

P30 Circulation relaxation parameter P30>1 (Use 1.0) 

BETAO Dcunping parameter 0<BETA)<0.25 (Use 0.15) 

FHALF: Switch to grid halving. 

=1. Halves the grid and need to provide Card 4a. 

=0. Does not use another grid (no need for Card 4a). 

FITMAX, FITMIN, CONV, PlO, P20, P30, BETAO 
Used for successively finer grids. 

Comments 

KTA5, PSIDEG, ALPHAS, TEMPF, NDTF, INSTFL, NGRD, LDTF 

KTA5: Forward flight speed, in knots 

PSIDEG: Azimuth location, «|» in deg (60*<4»<120®) 

ALPHAS; Shaft (or tip path plane) tilt in deg (+ aft) 

TEMPF: Outside air temperature in **F 

NDTF: DTF option 

=0, no DATAMAP connection 

^0, write DTF to connect to DATAMAP 

INSTFL; Instruction records option 

=0, no instruction records included for DTF 
^0, include instruction records for DTF. 

The last instruction record must be END and it must begin 
in column 1. 


80 



NGRD: Output grid option for DTF 

=0, use the default grid (appropriate for OLS comparison) 
^0, user will input the grid for DATAMAP use 

LDTF: File number for DTF. A positive value will write the DTF 

in internal format and a negative value will write the DTF 
in external format. 

=0, program resets it to the default value of 14 

CARD 7 ; Comments 

CARD 8 ; TIPSPD, RAD, RROOT, CR, CPOFFB, WAKPRT 
TIPSPD: R in ft/ sec 

RAD: Radius in inches 

RROOT: Root radius in inches 

CR: Reference chord for nondimensionalizing (in inches) 

CPOFFB : Controls printout of CP beyond blade trailing edge 0 : 

does not print beyond trailing edge 
#0: prints beyond T.E. (Use 1.0) 

WAKPRT: Controls print out of wake calculations (Use 0.0) 

CARD 9 : Comments 

CARD 10 : GFACTl, XVl, YVl, THVll, THV12 

CARD 11 : Comments 

CARD 12 : GFACT2, YV2, ZV2, THV21, THV22 

GFACTl and GFACT2 control the use of two infinite vortex line fila- 
ments in the flow field. 

XVl, YVl, THVll, and THV12 define the location and orientation of 
the first filament and similarly YV2, ZV2, THV21, and THV22 control 
the location and orientation of second filament in the flow field. 
It is recommended not to exercise these options. 

Use GFACTl =0.0 and GFACT2 =0.0 . 

CARD 13 : Comments 

CARD 14 : KPLOTS: Determines what grid stations along the blade you want the 

Cp, Mach Number, U, V, and W to be printed out. Typi- 
cally, the approximate blade root station is grid Sta #5 
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If NDTF=0 
If INSTFL 
CARD 15 ; 
CARD 16 ; 

CARD 17 : 
CARD 18: 


CARD 19 : 
CARD 20: 


If NGRD=0 
CARD 21 : 
CARD 22: 


and tip Sta is #21. You can print any desired station by 
writing the number of that station on this card. For the 
DATAMAP option, the grid stations included here must 
bracket the DATAMAP radial grid (FORMAT 2613). 

, omit cards 15 through 26. 

=0, omit cards 15 and 16. 

Comments 

DATAMAP instruction records. There can be »q) to 56 of these cards. 

The last one must be an END card (beginning in Column 1) (FORMAT 

18A4) . 

Comments 

GWHT, HCG, PSTAIC 

GWHT: Helicopter gross weight. Only used to label some DATAMAP 

plots. 

HCG: Helicopter Center of Gravity. Only used to label some 

DATAMAP plots. 

PSTAIC: Ambient air pressure. Used in CP to local Mach number 

calculations (Ib/in **2). 

Comments 

MODNUM, SHPNUM (FORMAT A8, 2X, A8) 

MODNUM: Helicopter Model Number (8 character field). Only used to 

Icdsel some DATAMAP plots. 

SHPNUM: Helicopter Ship Number (8 character field). Only used to 

lcd>el some DATAMAP plots. 

, omit cards 21 through 26. 

Comments 

NRADUS, NCHRD (FORMAT 213) 

The maximum allowed value of (2* NRADUS*NCHRD+1) is 510. The de- 
fault values are NRADUS=8 and NCHRD=20. 

NRADUS: Number of radial stations in the grid. If NRADUS=0, the 

radial stations are not redefined (the default OLS values 
are used). The maximiim value is 21. 
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NCHRD; Number of chord values in the grid. If NCHRD=0, the chord 
values are not redefined (the default values are used). 
The maximum value is 21. 

CARD 23; Comments 


If NRADUS 

=0, omit card 24. 

CARD 

24; 

DTFRAD 




DTFRAD; 

Radial values in fractions of the total radius (0 to 1). 
Input NRADUS values (7 per card) . The default values are 
(.40, .60, .75, .864, .910, .955, .970, .990). 

CARD 

25; 

Comments 


If NCHRD= 

0, omit card 26. 

CARD 

26; 

DTFCRD 

\ 



DTFCRD; 

Chord values in fractions of the total local chord (0 to 
1). Input NCHRD values (7 per card). The default values 
are (.005, .01, .02, .03, .05, .08, .15, .20, .25, .28, 
.32, .35, .37, .40, .45, .50, .55, .60, .70, .92). 

CARD 

_27; 

Comments 


CARD 

_28; 

FNC, PITCH, TAUP, XPO, BZ 



FNC: 

No. of spanwise stations used to define blade. 



PITCH: 

Blade pitch angle (cyclic collective) in deg + nose up. 



TAUP: 

Sweep of pitch axis w.r.t. ref. line, aft. 



XPO; 

Offset of pitch axis from ref line at Z=0 + toward T.E. 



BZ: 

Mesh spacing functions use of 1 gives uniform spacing. 0< 
BZ <1.0. (Use 0.389 for close spacing near tip.) 

CARD 

_29: 

Comments 


CARD 

30; 

ZZ, CHRD, 

THICK, XT, TWIST, XOFF, YOFF, FSEC (see Figure 50) 



ZZ: 

Span location in inches at which the following information 
is supplied. 



CHRD: 

Chord in inches at this ZZ. 



THICK: 

Linear scale factor. Multiplies with thickness ratio of 


airfoil defined in cards 30d. 
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XT: Location of twist axis in (x/CHD) from L.E. typical value 

is 0.25. 

TWIST; Geometric twist at this ZZ. 

XOFF: X=Location of L.E. w.r.t. ref line (+ toward T.E.) 

YOFF: Y=Location of L.E. w.r.t. ref line (+ upward). 

FSEC: =0 => No new airfoil defined. Uses the same airfoil as 

the one at the previous ZZ. 

=1 -> Reads a new airfoil at this ZZ. At grid locations 
between this zz and the previous ZZ, a cubic spline 
routine interpolates the airfoil geometry. 

If FSEC=1, provide cards 30a, b, c, and d. 

CARD 30a ; Consents 

30b: YSTM, FNU, FNL 


YSTM: 

1.0 

0.0 

for 

for 

symmetric airfoil (FNU=FNL) 
nonsymmetric airfoil 



FNU; 

No. 

of 

upper surface airfoil 

geometry 

def 

locations . 

FNL; 

No. 

of 

lower surface airfoil 

geometry 

def 

locations . 


30c ; Comments 
30d ; 1) XP, YP 

Provide FNU pairs of x/c and y/c coords of upper surface of the 
airfoil, one pair per card. 

If YSYH=1, then skip the remaining details of Card 30. 

If YSYM=0, then provide the following cards. 

2) CoRonent card 

3) XP, YP 

Provide FNL pairs of x/c and y/c coords of the lower surface of 
the airfoil, one pair per card. 

CARD 31 ; Repeat CARD 29 thru CARD 30 until FNC number of spanwise stations 
are defined. 


Sample input to R0T22 are provided in Tables 5 and 6 where the same blade 
geometry is provided in two different forms (see Figure 51). Of course, both 
of them will provide identical results. Sample input utilizing the DATAHAP 
interface is provided in Table 7. 
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TABLE 5. SAMPLE INPUT TO ARSA05 


TEST BLADE WITH NACA 0012 A/F&COMPLEX G!IOM«P • A^REF 


FNX 

FNY 

FNZ 

FCUNT 


8. 

12. 

0. 

F ITMAX 

FITMIN 

CONVO 

FHALF 

50. 

1 . 

0.00003 

1. 

40« 

!• 

0.00001 


FF .KNOTS 

PSl.OEG 

ALPHA. OEG 

AINFFPS 

100.0 

90.0 

-10.0 

1117.0 

TIPSP.FPS 

RAO. IN 

RROOT.IN 

REFCHORO 

700.0 

200.0 

50.0 

10.00 

K PLOTS 
lO 13 IS 

16 17 18 

19 20 21 


FNC 

PITCHOEG 

PA TILT 

PA OFF 

6* 

10.0 

2.0 

0.0 

STNvZZ 

CHORD 

TA/CHORO 

TWIST 

50.0 

10.0 

0.25 

5.0 

YSrN 

FNU 

fnl 


1*0 

44. 

44. 



rOFF 

0.0 


NACA0012 UPPER SURFACE 


0.0 

0.0 


0.00354 

0.01033 


0.00743 

0.01478 


0.01256 

0.01898 


0.01888 

0.02297 


0.02633 

0.02677 


0.03489 

0.03039 


0.04452 

0.03382 


0.05519 

0.03706 


0.06690 

0.04012 


0.07961 

0.04299 


0^09332 

A. 0AH6C» 


0.10800 

0.04612 


0.12364 

0.05037 


0.14024 

0.05241 


0.15777 

0.05421 


0.17623 

0.05578 


0.19562 

0.05712 


0.21591 

0.05820 


0.23710 

0.05904 


0.25919 

0.05962 


0.28216 

0.05994 


0.30601 

0.06001 


0.33073 

0.05961 


0.35632 

0.05935 


0.38276 

0.05863 


0.41006 

0.05764 


0.43821 

0.05639 


0.46719 

0.05489 


0.49702 

0.05313 


0.52768 

0.05111 


0.55916 

0.04884 


0.59147 

0.04633 


0.62460 

0. 04356 


0.65854 

0.04055 


0.69329 

0.03729 


0.72884 

0.03378 


0.76520 

0.03000 


0.80236 

0.02597 


0.84031 

0.02166 


0.87905 

0.01705 


0.91859 

0.01214 


0.95890 

0.00668 


1.00000 

0.00126 


ZZ 

CHORD 

XT 

IbO.O 

10.0 

0.25 

176.0 

O 

. 

o 

.4 

0.25 

179.0 

o 

• 

o 

0.25 

190.0 

e*o 

0.25 

200.0 

6*0 

0.25 


TWIST 

XOFF 

YOFF 

4.0 

-5.0 

0.0 

3.2 

-5.0 

O 

. 

o 

3.05 

-5.0 

o 

• 

o 

2.5 

1 

• 

O 

-1.0 

2.0 

-1.0 

0 
• 

(M 

1 


NEWSEC 

1 * 


NEWSEC 

0.0 

0*0 

0.0 

0*0 

0*0 
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TABLE 6 


SAMPLE INPUT TO ARSA05 


TEST BLADE WITH NACA 0012 A/FtCOMPLEX GEUM 


FNX 

FNY 

FN2 

FCONT 

60* 

6 * 

12* 

0* 

fitmax 

FITMIN 

CONVO 

FHALF 

so* 

!• 

0*00003 

1* 

40* 

1* 

0*00001 


FF* KNOTS 

PS 1* DEG 

ALPHA *OEG 

AINFFPS 

100*0 

90*0 

-10.0 

1117*0 

T APSP#FPS 

PAD* IN 

RRUUT .IN 

HEFCHORO 

700*0 

KPL0T5 

200*0 

50*0 

10*00 

lO 13 IS 

IS 17 18 

19 20 21 


FNC 

PITCHOEG 

PA TILT 

PA OFF 

6* 

10*0 

2*0 

5*0 

STN*ZZ 

CHORD 

TA/CHORO 

TWIST 

50*0 

10*0 

0*25 

5*0 

YSYM 

FNU 

FI4. 


1*0 

44* 

44* 


NACA0012 

0*0 

UPPER SURFACE 
0*0 



VOFF 

0*0 


0*00354 

0*01033 


0*00743 

0*01476 


0*01256 

0*01898 


0*01688 

0*02297 


0*02633 

0*026 77 


0*03489 

0*03039 


0*04452 

0*03382 


0*05519 

0*03706 


0*06690 

0*04012 


0*07961 

0*04299 


0*09332 

0*04566 


0*10800 

0*04812 


0*12364 

0*05037 


0*14024 

0*05241 


0*15777 

0*05421 


0*17623 

0*05578 


0*19562 

0*05712 


0*21591 

0*05820 


0*23710 

0*05904 


0*25919 

0*05962 


0*28216 

0*05994 


0*30601 

0*06001 


0*33073 

0*05981 


0*35632 

0*05935 


0*38276 

0*05863 


0*41006 

0*05764 


0*43621 

0*05639 


0*46719 

0*05489 


0*49702 

0*05313 


0*52768 

0*05111 


0*55916 

0*04884 


0*59147 

0*04633 


0*62460 

0*04358 


0*65854 

0* 04055 


0*69329 

0*03729 


0*72884 

0*03378 


0*76520 

0*03000 


0*80236 

0*02597 


0*64031 

0*02168 


0*67905 

0*01705 


0*91859 

0*01214 


0*95690 

0*00688 


1*00000 

0*00126 


ZZ 

CHORD 

XT 

150*0 

10*0 

0*25 

176*0 

o 

• 

o 

PR 

0*25 

179.0 

o 

• 

o 

PR 

0*25 

190*0 

8*0 

0*25 

200*0 

6*0 

0*25 


IWIST 

XOFF 

YOFF 

4*0 

0*0 

0*0 

3*2 

o 

• 

o 

0*0 

3*05 

o 

• 

o 

O 

• 

o 

2*5 

o 

• 

CM 

-1*0 

2*0 

4*0 

-2*0 


NEWSEC 

1 . 


NEWSEC 

0»0 

0*0 

0«0 

0*0 

OmO 
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Figure 51. Blade geometry definition. 


TABLE 7. SAMPLE INPUT TO ARSA05 


OLS AH-IG DATA CREATION WITH ^ 

ROT22 FOR 

DATAMAP 

INPUT 


FNX 

FNY 

FNZ 

PLTCP 

PLTSEC 

PLTSRF 

FCONT 

30. 

4. 

8. 

0. 

1 . 

1 . 

0. 

FITMAX 

FITMIN 

CONV0 

P10 

P20 

P30 

BETA0 

10. 

1. 

0.00001 

1.6 

1 .0 

1 .0 

0.15 

20. 

1. 

0.00001 

1.725 

1 .0 

1 .0 

0.15 

30. 

1. 

0.00001 

1.80 

1 .0 

1 .0 

0.15 

KTAS 

PSI,DEG 

ALPHA, DEG 

TEMF 

NDTF 

INSTFL 

NGRD 

150.2 

90.0 

-4.329 

68.9 

1 . 

1 .0 

0.0 

TIPSPD 

RADIUS 

RROOT 

CR 

CPOFFB 

WAKPRT 


694.0 

264.0 

40.0 

28.63 

1 .00 

0.0 


GFACTl 

XVI 

YVl 

THVll 

THV12 



0.0 

0.25 

0.25 

10.0 

5.0 



GFACT2 

XV2 

YV2 

THV21 

THV22 



0.0 

-0.5 

12.1 

0.0 

0.0 



KPLOTS 







11 12 13 

14 15 16 17 18 19 20 

21 22 23 

IS) 

IS) 

cn 

ro 

27 28 29 

30 31 32 

DATAMAP 

INSTRUCTION 

RECORDS BEGIN 




REP ROT2872 

PASSWRD TOM 




SPACE 

2000 

USERNAM K 

.R.SHENOY 




ITEMS 

- 






NOFILT 

CAL IB 

RATE 

100 





ALL ITEMS 
COUNTERS 

RECORD 10 ABSOLUT 0 

ALL COUNTERS 
END 


GW 

HCG 

PSTAT 




8500. 

0.0 

11.98 




MODEL NO 

SHIP NO. 





OLSAHl-G 

REC 2872 





FNC 

COLP 

TAUP 

XP0 

BZ 


8. 

12.185 

0.0 

6.88 

0.389 


STN,2Z 

CHRD 

THICK 

XT 

TWIST XOFF 

YOFF 

40.0 

28.63 

1 .0 

0.2403 

-2.25 0.0 

0.0 

YSYM 

FNU 

FNL 




1 . 

48. 

48. 




XP 

YP 

: STA 

40(REF CHORD 

IS 28.3 IN.) 


0.0 

0.0 






0.00453 0.01163 
0.00862 0.01569 
0.01250 0.01857 
0.01628 0.02087 
0.02366 0.02453 
0.03092 0.02743 
0.03810 0.02986 
0.04524 0.03195 
0.05590 0.03463 
0.06652 0.03689 
0.08064 0.03939 
0.09121 0.04098 
0.10177 0.04235 
0.11231 0.04353 
0.12285 0.04455 
0.13339 0.04543 
0.15093 0.04660 
0.17548 0.04773 
0.20001 0.04836 


FPRINT 
-1 . 
FHALF 
1 .0 
1 .0 

LDTF 

14. 


FSEC 

1 . 
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TABLE 1 . (Concluded) 


0.21402 0.04851 
0.22102 0.04855 
0.22452 0.04855 
0.23503 0.04851 
0.25253 0.04832 
0.28053 0.04767 
0.30152 0.04696 
0.32950 0.04572 
0.35048 0.04461 
0.38544 0.04244 
0.42036 0.04015 
0.45528 0.03785 
0.49020 0.03556 
0.52512 0.03326 
0.56003 0.03097 
0.59495 0.02867 
0.62987 0.02638 
0.66479 0.02408 
0.69970 0.02179 
0.73462 0.01949 
0.76954 0.01720 
0.80446 0.01490 
0.83938 0.01261 
0.87430 0.01031 
0.90921 0.00801 
0.94413 0.00572 
0.97905 0.00342 
1.00000 0.00205 


STN,2Z 

80.0 

CHRD 

28.63 

THICK 

1.0 

XT 

0.2403 

132.0 

28.63 

1.0 

0.2403 

198.0 

28.63 

1.0 

0.2403 

224.0 

28.63 

1 .0 

0.2403 

238.0 

28.63 

1.0 

0.2403 

251 .0 

28.63 

1.0 

0.2403 

264.0 

28.63 

1 .0 

0.2403 


TWIST 

XOFF 

YOFF 

FSEC 

-4.11 

0.0 

0.0 

0. 

SEC 

-6.4 

0.0 

0.0 

0. 

SEC 

-9.14 

0.0 

0.0 

0. 

SEC 

-10.2 

0.0 

0.0 

0. 

NGOP 

-10.74 

0.0 

0.0 

0. 

NGOP 

-11.26 

0.0 

0.0 

0. 

SEC 

-11.78 

0.0 

0.0 

0. 



OUTPUT DESCRIPTION 


The first section of the output describes the geometry of the blade. At each 
section Z, the spanwise location is printed out in both ZZ/CREF and in inches. 
The XSING and YSING correspond to the singular point near the leading edge of 
the blade about which a transformation is made. For practical purposes, the 
user can assume that the singular point is at the L.E. After all the sections 
are printed, ref. chord, sweep at the root of the blade (SWEEPl), sweep at the 
tip (SWEEP2), and the sweep of the singular line at infinity (SWEEP is set= 
SWEEP2/2) are printed. 

Similarly, the dihedral angles are printed. The grid information and other 
transformation details follow. Finally, the section characteristics are 
printed out for the rough grid. The finer grid section characteristics are 
those of interest and they include freestream Mach No. at the section, azimuth 
location, shaft tilt angle, spanwise location (r/R), section c., , c^ and 
^ml/4c' ^pcr chord. Further, the x locations at the given section in 

inches, y/c at that x, local surface Mach number, and U, V, W components of 

surface velocities are printed out. A graphical representation of the 
values is provided along a distorted X scale. The location of the leading 
edge is also indicated. The information above the leading edge location cor- 
responds to the lower surface 2 md below, the leading edge location corresponds 
to the upper surface. Following this, a listing of sectional characteristics 
and their blade characteristics is printed. Finally, the DATAMAP file crea- 
tion program details are printed. 

User Tips 

1. If there are sharp changes (kinks) in the planform, the interpolating 
cubic spline will exhibit instabilities unless the following steps are 
taken: 

i) Do not specify any geometry information at the kink. 

pfscmm 9^ mi umbo 
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ii) Define geometry on either side of the kink (r/R=±0 .005) . However, 
if the blade has just a single, straight sweep, one need not be 
concerned about the location of the point after the kink (near the 
tip). 

iii) Define at least two points after each kink. 

For a swept blade with smooth plcmform, define a number of stations to 
represent the blade as accurately as possible. However, if the airfoil 
definition changes from section to section, one may be limited by the 
number of stations at which the airfoil information is available. One 
may try to use the "thick" pareuneter to define transition sections. 

The information on the blade surface is calculated at the computational 
grid points located evenly in the transformed plane. Consequently, the 
results are printed at nonuniform intervals in the physical plane. 

The program uses a left-hand sign convention. As a result, x is in the 
chordwise direction, positive toward the TE; y is perpendicular to the 
chord and is positive upwards; z is in the spanwise direction, positive 
toward the tip. 
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